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Braneworld gravity is a model that endows physical space with an extra dimension. In the type II 
Randall- Sundrum braneworld gravity model, the extra dimension modifies the spacetime geometry 
around black holes, and changes predictions for the formation and survival of primordial black holes. 
We develop a comprehensive analytical formalism for far-field black hole lensing in this model, using 
invariant quantities to compute all the geometric optics lensing observables: bending angle, image 
position, magnification, centroid, and time delay. We then make the first analysis of wave optics 
in braneworld lensing, working in the semi-classical limit. Through quantitative examples we show 
that wave optics offers the only realistic way to observe braneworld effects in black hole lensing. We 
point out that if primordial braneworld black holes exist, have mass M», and contribute a fraction 
fhh of the dark matter, then roughly ~ 3 x 10 5 x /bh(M./10 -18 M Q ) _1 of them lie within our 
Solar System. These objects, which we call "attolenses," would produce interference fringes in the 
energy spectra of gamma-ray bursts at energies E ~ 100 (A/./10" 18 M©) -1 MeV (which will soon 
be accessible with the GLAST satellite). Primordial braneworld black holes spread throughout the 
universe could produce similar interference effects. If they contribute a fraction fi. of the total energy 
density, the probability that gamma-ray bursts are "attolensed" is at least ~0.1f2.. If observed, 
attolensing interference fringes would yield a simple upper limit on M. . Detection of a primordial 
black hole with M, < 10 -19 M@ would challenge general relativity and favor the braneworld model. 
Further work on lensing tests of braneworld gravity must proceed into the physical optics regime, 
which awaits a description of the full spacetime geometry around braneworld black holes. 

Keywords: gravitational lensing, gravity theories, extra dimension 



I. INTRODUCTION 

Gravitational lensing has emerged as a powerful and 
far-reaching tool in astrophysics and cosmology Q, 0, Q . 
In this series we are showing how lensing can also be 
employed to test theories of gravity. Papers I |4j and 
II [5j developed an analytical formalism for identifying 
the lensing signatures of gravity models that fall within 
the post-post-Newtonian (PPN) framework, even prob- 
ing out to third order in such models. These studies un- 
covered some surprising universal relations among lensing 
observables that helped us make specific predictions that 
are testable with current or near-future instrumentation. 

In this paper we examine a gravity model that lies 
outside the PPN framework, namely, type II Randall- 
Sundrum braneworld gravity [jj. According to this 
model, familiar 4-dimensional spacetime is actually a 
submanifold (a "brane") in a 5-dimensional spacetime 
(the "bulk"), with the extra dimension characterized by a 
curvature radius £ which could be as large as ^0.2 mm |7|. 
One intriguing prediction of the model is that braneworld 
black holes might be produced at energies as low as ~1 



TeV, which could lead to observable Hawking radiation 
in the forthcoming Large Hadron Collider Q or create 
specific signatures in cosmic ray showers 0. 

Another important prediction is that braneworld black 
holes produced in the early universe might survive to 
the present day. Primordial black holes are predicted to 
have formed from density fluctuations in the very early 
universe, with a mass spectrum that increases rapidly 
towards low masses [id. Il II . In general relativity, black 
holes smaller than ~ 1CP 19 Mq would have evaporated 
by now through Hawking radiation |12| . Compared with 
their GR counterparts, however, braneworld black holes 
evaporate more slowly 01 an d ma Y have accreted more 
efficiently in the early universe [Til Il5j . Together, these 
effects may allow primordial braneworld black holes as 
small as M. ~ 1 kg ~ 10 -30 M Q to survive to today 
[3G3. The implication is that primordial braneworld 
black holes may contribute some fraction of the unknown 
dark matter. Gravitational lensing offers a crucial test 
of this hypothesis if we can identify appropriate lensing 
scenarios that carry a clear imprint of braneworld gravity. 
That is our goal. 
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An exact metric for the spacetime geometry induced 
by a braneworld black hole is still unknown. In the 
far-field or weak-deflection regime, it is well established 
that braneworld black holes are described by the Garriga- 
Tanaka metric [13, Il8l Il9| . In the near-field or strong- 
deflection regime, various metrics are believed to approx- 
imate the true spacetime geometry (e.g., [2(j). These 
different metrics have been used to explore gravitational 
lensing by braneworld black holes (sec the review by Ma- 
iumdar & Mukherjee 0). For example, Kar & Sinha 
|2ll | computed the light bending angle for the Garriga- 
Tanaka and tidal Reissner-Nordstrom metrics. Majum- 
dar & Mukherjee determined the light bending an- 
gle, image position, and magnification for the weak- 
deflection regime of the Myers-Perry metric. Eiroa |2^| 
and Whisker [24| studied the bending angle, image posi- 
tion, and magnification for the strong-deflection regime 
of the Myers-Perry and tidal Reissner-Nordstrom met- 
rics, respectively. Still other metrics may prove use- 
ful for studying the approximate lensing properties of 
braneworld black holes (cf. [MEl). 

We present a thorough study of lensing in braneworld 
gravity, including a realistic assessment of prospects for 
observing braneworld effects in astrophysical lensing sce- 
narios. We focus on weak-deflection lensing for two 
reasons. First, the images that appear in the strong- 
deflection limit carry important near-horizon effects but 
are exceedingly difficult to observe [2{||2(j. Second, we 
shall argue that wave optics will play a crucial role in 
lensing tests of braneworld gravity, and wave optics ob- 
servables are dominated by the two images that appear 
in the weak-deflection regime. Since the Garriga-Tanaka 
metric correctly describes the spacetime geometry in the 
far- field regime of a braneworld black hole, we develop 
the full analytical infrastructure for weak-deflection lens- 
ing in this metric. 

We rederive the light bending angle in the Garriga- 
Tanaka metric, but express it for the first time in terms 
of invariant quantities (Section [nj. We go beyond the 
bending angle and use invariant quantities to compute 
the observable properties of the lensed images: positions, 
magnifications, and time delays (Sections lllllll V|) . We 
then consider for the first time wave optics effects in 
braneworld black hole lensing (Section 0) . Finally, we 
examine a variety of applications of braneworld lensing. 
We show that traditional astrophysical lensing scenar- 
ios will be unable to measure braneworld effects in the 
foreseeable future (Section IVI|) . However, an application 
of wave optics that we call "attolensing" provides excit- 
ing opportunities for observing braneworld effects (Sec- 
tion [^Rj. In particular, we point out that if primordial 
braneworld black holes exist and contribute to the dark 
matter, they must exist not only throughout the universe 
but also within our Solar System. Attolensing will pro- 
vide a crucial test of this prediction of the braneworld 
model. 
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FIG. 1: Schematic diagram of the lensing geometry. Stan- 
dard quantities are defined as follows: B is the angular po- 
sition of the unlensed source; 1? is the angular position of an 
image; a is the bending angle; and dh, ds, and dhs are angu- 
lar diameter distances between the observer, lens, and source. 
The impact parameter b is an invariant of the light ray and is 
related to the angular image position by # = sin~ 1 (&/dz,). 

II. METRIC AND LIGHT BENDING ANGLE 

We begin by stating our basic assumptions (cf. 0,0)- 
Consider a gravitational lens with mass M, that is com- 
pact, static, and spherically symmetric, with an asymp- 
totically flat spacetime geometry sufficiently far from the 
lens The spacetime is vacuum outside the lens and 
flat in the absence of the lens. We adopt the standard 
lensing scenario shown in Figure ^ with the observer 
and source lying in the asymptotically flat regime of the 
spacetime. 

The light ray's distance of closest approach ro and im- 
pact parameter b are both assumed to lie well outside 
the lens's gravitational radius m. = GM,/c 2 . The light 
bending angle is assumed to have the following form at 
lowest order in m.: 

where A\ and B2 are independent of m. and b. Since 
b and m. are invariants of the light ray, this expression 
for the bending angle is independent of coordinates. If 
B2 ^ 0, the bending angle cannot be written as a 
series in the single parameter m,/6, which places this 
model outside the PPN framework studied in Papers I 
and II. We now show that braneworld black holes do 
yield bending angles of this form. 

A. Isotropic coordinates 

The Garriga-Tanaka metric is often written as follows 
in isotropic coordinates: 
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This form is valid only in the limit 



m. 



— < 1, < 1, 



(3) 



and the exact metric describing the spacetime geometry 
around braneworld black holes is not yet known. Wc 
shall verify a posteriori that our lensing solutions satisfy 

The metric © has terms of the form ±2m,/r as in the 
weak-deflection limit of general relativity. Braneworld 
effects enter via the I terms. (When I — we recover the 
standard far-field black hole metric of general relativity.) 
Notice that the I terms factor as (m,/f)(£ 2 /f 2 ), so we 
can adopt the approach of taking Taylor series expansions 
in m./f to obtain the appropriate weak-deflection limit 
including the braneworld terms. Later we shall consider 
when it is appropriate to take series expansions in Ijf as 
well. 

We can think of the metric © more generally as having 
the form 

ds 2 = -A(f) alt 2 + B(f) df 2 + C(f ) f 2 dn 2 , (4) 
where the metric functions A(f) and B(f) — C(f) are 



readily identified. For such a metric, the distance of clos- 
est approach fo is related to the impact parameter b by 
(cf. eq. 12 of Paper I) 



b = r 



C(fo) 



-,1/2 



MM 
i 



m 



Inverting this relation yields 



ro 



b 2 



m. 

b 







m. 



(5) 



(6) 



For a metric of the form the light bending angle 
can be written as 



d(r ) 



1 



A B 



-,1/2 



Ci/V-AC/r 2 



dr 



7T. (7) 



Plugging in the metric functions, and temporarily replac- 
ing b with fo using 0, we can write a series expansion 
for the integrand, 



ro 



r (r 2 



})l/2 



1 + 



2f 2 f 2 l + e 2 (f 2 + ff + fo) m « 



rr 2 (r + r ) 



O 



ro 



m. 



dr 



(8) 



Carrying out the integration yields the deflection angle 
in terms of the isotropic coordinate distance of closest 
approach, 



(9) 



This agrees with the bending angle found by Kar & Sinha 
21]. However, as emphasized in Paper I, expressions like 
Q are coordinate dependent and should be re-expressed 
in invariant form. Using © to rewrite the distance of 
closest approach fo in terms of the invariant impact pa- 
rameter b yields 



a(b) = 4 



m. 



r 2 
b 2 



o 



(10) 



At this order of approximation the form of the bending 
angle is the same for ro and 6; still, it is important to 
use the invariant expression. Notice that IjlQI) has the 
form assumed in QJ, with = 4 and B 2 = £ 2 , which 
shows that braneworld black hole lensing lies outside the 
standard PPN framework. 



B. Standard coordinates 

As a consistency check of (|1U|I . and to connect with 
our formalism in Papers I and II, we rederive the bend- 
ing angle starting from standard coordinates (or the area 
gauge 20] ). In this case we write the metric in the form 

ds 2 = -A(r)dt 2 + B{r)dr 2 +r 2 dQ 2 . (11) 

Comparing the dVt 2 terms in and (|11|) shows that the 
isotropic and standard radial coordinates are related by 









1 + 










3f 2 ) 



m. 



O 



m 



Inverting this relation yields 

I 2 



1-1 



3r 2 



m. 



(12) 



(13) 



Returning to the metric (J2J and changing radial coordi- 
nates yields a metric of the form (flip with metric func- 
tions 



1-21 



A(r) 

B(r) = 1 + 211 



6r z I r \ r 



(14) 
(15) 
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When the metric is written in the form the The light bending angle can be written as 

standard-coordinate distance of closest approach r and 
invariant impact parameter b are related by 



b = 



ro 



ro 



I, w\ 













ro 



m. 



2£ 2 \ m. m /m.^ 



(16) 



a{r ) = 2 



1 



AB 



-,1/2 



1/b 2 - A/r 2 



dr - ir. (18) 



Plugging in the metric functions, and temporarily replac- 
ing b with ro using l|17(l . we can write a series expansion 
(17) for the integrand, 

I 



a(r ) = 2 



r (r 2 — r 2 ) 1 / 2 



3r 2 r 2 (r 2 + rr + r%) + ^ 2 (2r 4 + 2r 3 r + 2r 2 r 2 , + 3rrj] + 3r^) m. 



01^ 

ro 



dr 



3r 3 r 2 (r + r ) 



ro 



(19) 



Carrying out the integration yields the deflection angle 
in terms of the coordinate distance of closest approach, 



ro 



ro 



(20) 



Using l|17|) to rewrite the distance of closest approach r 
in terms of the impact parameter b yields 



m, 



(21) 



This result agrees with (tTU|) , showing that we arrive at the 
same desired invariant bending angle expression starting 
from isotropic and standard coordinates. 



III. IMAGE POSITIONS, MAGNIFICATIONS, 
AND CENTROID 



We now go beyond the bending angle to determine 
observable quantities in braneworld black hole lcnsing. 
We examine the image positions and magnifications in 
this section, and defer the time delay to Section lTVl (This 
analysis parallels Section IV in Paper I.) We begin with 
the general lens equation (cf. Figure ^1, 



tan£ = tani? - D (tantf + tan(a - #)), 



(22) 



where B is the angular position of the source, § — 
sin _1 (6/dL) is the angular position of the image, and 
D = d^s/ds- We shall see that the lens equation yields 
two images in the far-field (weak-deflection) regime, one 
on the same side of the lens as the source and the other 
on the opposite side. Following the convention of Papers 
I and II, angles describing image positions are taken to 



be positive. This forces the source's angular position to 
have different signs: B is positive when the image is on 
the same side of the lens as the source (as depicted in 
Figure while B is negative when the image is on the 
opposite of the lens from the source. 

We now seek an appropriate series expansion of the lens 
equation in the weak-deflection limit. First, we change 
variables in light of the fact that lensing quantities natu- 
rally scale with the weak-deflection angular Einstein ring 
radius, 



4GM.d LS 

c 2 dLds 



Specifically, we define: 



B_ 



AD ' 



(23) 



(24) 



where i9. = tan _1 (m./<ii). In other words, the quanti- 
ties (3 and 9 are the scaled angular positions of the source 
and image, respectively. The quantity e m represents the 
angle subtended by the gravitational radius normalized 
by the angular Einstein radius, and it replaces m,/6 as 
our expansion parameter. We also define an angle asso- 
ciated with the braneworld length scale £, 



•dt = tan 



and a scaled version of this angle 



(25) 



(26) 



When it becomes appropriate to consider series expan- 
sions in Ijb (below), we shall actually use ei as the ex- 
pansion parameter. The conditions for validity of the 
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Garriga-Tanaka metric are equivalent to the conditions 
e. m <C 1 and £ e <C 1. 

With these substitutions, and the bending angle from 
(|2ip. the lens equation becomes 







-p + e 



i 



-2 \ 1 



1 



E m + O (S m y 



(27) 



To determine braneworld effects on the image positions, 
we then postulate that the position can be expanded in 
the form 



V = V + ViEi + V 2 e 
Plugging this into l|27[l yields 
1 



0(e e ) 



= 



-f 3 + ^-Q- Q 







) + (- 









(28) 



(29) 



(1 + gg)flofl 2 - 1 - 



9 2 

— Ee 



This is the desired series expansion of the lens equation. 

The zeroth-order term in l|29l) is just the usual lens 
equation for the weak-deflection limit of general relativ- 
ity, whose solution is 



■0 



(30) 



We neglect the negative solution because by conven- 
tion angles describing image positions are positive. The 
positive-parity image 6q , which lies on the same side of 
the lens as the source, is found by using (3 > 0. The 
negative-parity image 6q is then found by using [3 < 0. 
Explicitly, 



(31) 



The first-order term in l|29(l can be satisfied only if 
9\ = 0. In other words, there is no braneworld correction 
to the lensed image positions at first order in e^ ~ £/b. 
This is not surprising, since braneworld effects enter the 
metric at order (£/r) 2 . 

The second-order term in (|29|) is satisfied if 



(32) 



Thus, the full expression for the image position in 
braneworld gravity is 



= c 



9o(l + #o 2 ) 



(33) 



Rewriting 9q in terms of (3 using (|31|l . we can express the 
image positions in terms of the source position as 



if 2 + (3 2 



(34) 



2 \ yiT^ 



O (e^ 



The coefficient of e\ is positive for all value of (3, which 
means that braneworld effects push both the positive- 
and negative-parity images farther away from the lens 
(relative to the results for general relativity). 

By spherical symmetry, the signed magnification /i of 
a lensed image at angular position i3 is 



sinB(tf) dB(6) 



smi3 



dO 



(35) 



After changing to our scaled variables, we first make a 
Taylor series expansion in s m : 



M = 



0(e r . 



(e i + e 2 + iE 2 W-8 2 ~E 2 ) 

Now using (|33|l for the image position, we find 



(36) 



01 



2 6^ 



9^-1 



9o 4 -l)(^o 2 + l) 3 



+ 0(EiY. (37) 



The absolute magnifications in terms of the source posi- 
tion are given by 



1 



/ 2 + /3 2 



± 1 



(38) 



2 Vl/?| v / 4+^ , 

" |/3|(4 + /32)3/ 2 £ ' + 0(£f)3 ' 

where fjT < 0. Observe that at lowest order the 
braneworld magnifications have /i + + fj," = 1. This is 
identical to the lowest-order universal magnification re- 
lation in eq. (36) of Paper II. It is interesting to see that 
braneworld gravity obeys the relation originally derived 
for PPN models, even though it lies outside the PPN 
framework. 

In cases where the two images cannot be separately re- 
solved (such as microlensing 0,00], H2), the useful quan- 
tities are the total magnification and the magnification- 
weighted centroid position. The total magnification can 
be written as 



Mtot 



1 



\(3\^/JTW 1/31(4 + /? 2 ) 3 / 2 



(39) 
e\ + 0(e £ ) 3 . 



The magnification-weighted centroid position can be 
written as 



1/1 1 



0-\fM- 



'I (3 + /3 2 ) 
2 + /3 2 



(40) 



(2 + (3 2 ) 2 



0{et)- 



We see that braneworld effects decrease the total mag- 
nification and push the centroid farther from the lens 
(compared with the results for general relativity). 
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IV. TIME DELAY 

We now derive the lensing time delay (in parallel with 
Section V of Paper I) . We first focus on a spacetime that 
is static and asymptotically flat, and discuss the general- 
ization to a curved universe cosmology at the end of this 
section. 

Let i? src and R bs be the radial coordinates of the 
source and observer, respectively. From geometry rel- 
ative to the flat metric of the distant observer (who is 
assumed to be at rest in the natural coordinates of the 
metric eq. 1111) , we can work out (see Figure ^| 



Roh 



R s 



l LS 



4 tan 2 B) 1 ' 2 . 



(41) 



The radial distances are very nearly the same as angular 
diameter distances since the source and observer are in 
the asymptotically flat region of the spacetime. In other 
words, the distortions in distances near the black hole 
are assumed to have little impact on the total flat met- 



ric distance from the compact body to the observer or 
source. 

In the absence of the lens the light ray would travel 
along a linear path from the source to the observer with 
length ds/ cosB. The time delay r is the difference be- 
tween the light travel time for the actual ray, and the 
travel time for the straight line the ray would have taken 
in the absence of the lens. This can be written as 



cT = T(R sl . c )+T(R ohs 



with (see eq. 91 of Paper I) 



ds 
cosB 



(42) 



T(R) = - b 



1 / A(r)B{r) 
~A{r) V l/b 2 - A(r)/r< 



dr . (43) 



Using the metric functions l|14|l and (|15fl . we find for 
braneworld gravity 



T{R) = 



3r 2 r (2r + 3r ) + ^ 2 (2r 2 + 7rr a + 7rg) m. 



3r 3 ra(r + r ) 



R 2 



2 In 



R+ 



ro 



+ ^R-r (l + i 



9R + 7r 
3Rr 2 



O 



dr, 
m. 



(44) 



To obtain a coordinate-invariant expression, we could re- 
place ro with b using (fTTj) . We would then want to take 
a series expansion in b/ R as well as m,/&, because they 
are of the same order (see Section V of Paper I) . 

It is simpler to proceed directly to the time delay ex- 
pressed in terms of our scaled angular variables. We com- 
pute T(R src ) and T(i? b s ) using the radii from (|4*T)l . We 
change to angular variables using b = dz, smi9, and then 
reintroduce the scaled angular variables 9 and P defined 
in J23J). We work to lowest order in e m , and then take a 
Taylor series in the dimensionless braneworld parameter 
eg. The result is 



TE 



1 + (3 2 - el - In 



d L e, 



)2 Q2 
V E 



4d L s 

where the natural time scale is 



TE 



cd LS 



7? 2 



m. 



2# + °^ 3 ' 
(45) 



(46) 



Note that we have used l|33|l for the image position to 
obtain an expression written in terms of 6q, the image 
position in the weak-deflection limit of general relativity. 

More interesting than the individual time delays is the 
differential delay At = t~ — r + between the negative- 



and positive-parity images: 



At 

TE 



- (Oof 



ln^ 



vr 2 



+ 0(e e )\ (47) 



-1/31^4+^ 



In 



+ if 1/3174+^ + 0(e ( y 



(48) 



The first expression is written in terms of the image po- 
sitions 9^, while the second is written in terms of the 
source position. In each case, the order unity term re- 
covers the familiar time delay in the weak-deflection limit 
of general relativity. To simplify the notation below, we 
write the differential time delay as 



At 

— = A + A! + A e] + O {e t f 

TE 



(49) 



where 



Ao = i|/3|v / 4+^ 2 ", 



Ai = In 



(50) 
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It is straightforward to place this analysis in a back- 
ground universe that is curved and expanding. We are 
working in the limit that the impact parameter is small 
compared to the distances between the observer, lens, 
and source. Thus, the light path is determined by the 
background cosmology for all but the tiny fraction of the 
path when it is near the lens. (See [1 @ for further 
discussion.) All of the previous analysis holds if we in- 
terpret d.L, ds, and cIls as cosmological angular diameter 
distances, and we modify the natural time scale from i|46|) 
to 



te = (1 + z L ] 



cd LS 



4(1 



m. 

zl) — , 

c 



(51) 



where Zl is the cosmological redshift of the lens. 



V. SEMI-CLASSICAL INTERFERENCE OPTICS 

The foregoing analysis applies in the limit of geometric 
optics. We now examine how the wave nature of light 
can affect gravitational lensing observables. Wave optics 
have been studied for standard weak-deflection lensing in 
gen eral relativity (see the review by Nakamura & Deguchi 
[27j . and Sections 4.7 and Chapter 7 of Schneider et al. 
[lj), but to our knowledge have not been treated before 
in the braneworld lensing literature. 

In analogy with the Young double-slit experiment, 
we may consider that light waves from the positive- 
and negative-parity images interact to produce an in- 
terference pattern at the observer. To quantify differ- 
ent regimes, suppose a pointlike light source emits a 
monochromatic, spherical light wave of period T = 2n/u 
or wavelength A = cT, which is lensed by a braneworld 
black hole. Geometric optics apply in the limit T — > 0. 
If the period is finite but small compared with the time 
delay between the imag es ( T <C At), then we are in the 
semi- classical limit (cf. 27, 38]). This is the regime we 
shall investigate. Larger periods T > At lead to the 
physical optics limit. In this regime, all regions of the 
lens plane contribute to interference effects. The analysis 
therefore requires knowledge of the time delay function 
across the entire lens plane, which in turn requires knowl- 
edge of the metric at all such positions. However, the 
full metric for braneworld black holes is not yet known, 
which means that we cannot give a complete wave optics 
treatment of lensing in braneworld gravity at this time. 

There is a well-established connection between the ge- 
ometric optics and semi-classical limits. If geometric op- 
tics predict a set of images with magnifications fa and 
time delays Ar^ , then the total magnification in the 
semi-classical limit is given as follows (e.g., Sections 4.7 
and 7.1 of 1], and Section 2 of H3): 



j i<3 



(52) 



where = 0, 1, 2 depending on whether the fc-th image 
is a minimum, saddle, or maximum, respectively. The 
first sum gives the total magnification in the geometric 
optics limit, while the second captures the corrections 
due to semi-classical interference between the images. 
There is a phase factor of wAry from the lensing time 
delay, and an additional factor of (nj — 7ij)7r/2 from dif- 
ferences in the phases of the images upon exiting the lens 
plane. 

Lensing by a black hole produces a positive and neg- 
ative image pair, so we have ni2 = 1 — = 1 and 
Ari2 = r~ — t + . The total magnification in the semi- 
classical limit can then be written as 



M = \[i + \ + \fjT\ + 2|/iV~| 1/2 sin(wAr 



(53) 



The sine term creates a series of bright and dark fringes 
in the magnification as a function of energy or wave- 
length. Bright fringes (maxima in M) occur when the 
phase difference is wAr = 2jir + tt/2 for integer j, which 
correspond to photon energies 



E Y - (j + V4) ^ 



J = 0,1, 2,. 



(54) 



Dark fringes (minima in M) occur when the phase dif- 
ference is luAt = (2j + l)ir + tt/2, or energies 



E 



dk 



(i + 3/4)^ 



3 = 0,1,2, 



(55) 



Note that the energy spacing between adjacent bright 
fringes, or adjacent dark fringes, is always AE = h/Ar, 
independent of j. 

For a braneworld black hole, the magnification sum 
1 + | is given in (|39|l , and the magnification product 

is 



1 1/2 



1 



2 + f3 2 



\/3\y/I+W 1/31(4 + /3 2 ) 3/ 



2 e 2 e +0(e e ) 3 . 

(56) 

Consequently, the semi-classical total magnification is 
given as follows in terms of the source position: 



M 



2 + 1 + 2 sin(w At) 
4 + /3 2 



(57) 



2(2 + /? 2 )sin( W Ar) 2 3 
— e e + (J {ee) 



(4 + /3 2 ) 2 



Using the time delay from l|49|l . we find that the bright 
fringes have magnification 



M hl = 



1 - 



4 + /3- 



■4 + o(e e y 



(58) 



and are located at energies 

br C? + l/4)*> 



E 



(A + Ai)t e 



1 



A 



A + Ai 



(59) 
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while the dark fringes have magnification 



M dk = 



i + 



V4+^ 
and are located at energies 
0' + 3/4)A 



4 + /3 : 



2 c e 



E dk 

j (A + Ai)r B 



1 - 



A 



Ao + Ai 



0(e e ) 



4 + 0(eiY 



(60) 



(61) 

The absolute peak-to-trough distance between bright and 
dark fringes is 



M hx_ M dk = 



1 - 



4 + /3 



2 e l 



0(e e ) 



while the fractional difference is 



(62) 



M hl - M 



<ik 



M 



dk 



2 + /3 2 



1 - 



P 2 



2 + /3 



2 °t 



o(e t y 



We see that braneworld effects (from the extra dimension 
of space) shift the interference fringes to lower energies, 
and contract the energy spacing E^ k — Ej. Braneworld 
effects also reduce the magnification of the bright fringes 
and increase the magnification of the dark fringes, which 
reduces the peak-to-trough distance in the fringe pattern. 



VI. LARGE BRANEWORLD BLACK HOLES IN 
TRADITIONAL LENSING SCENARIOS 

We have seen that all the braneworld corrections to 
weak-deflection black hole lensing scale with 



tan _1 (£/dz 



0* 



AGM. d L d LS 4m. d L d LS 

(64) 

where we used tan 1 (l/di J ) « lj&L and the definition of 
•&E from H23|) . Notice that eg is given by the ratio of the 
braneworld scale I to the geometric mean of the gravita- 
tional scale m. and the astrophysical distance scale. 

There are two ways to maximize braneworld effects. 
One is to consider a small black hole mass, since e\ oc 
m^ 1 . The other possibility is to make either dj, or d^s 
small. Notice that if d^ <C d$ then d^s ~ ds and the 
combination of distances reduces to 1/di,. Alternatively, 
if dLS *C ds then dL ~ ds and the combination of dis- 
tances reduces to 1/<1ls- We can combine these two pos- 
sibilities and write 

1)2 

d min = mm(d L , d L s) < d s ■ (65) 



4 m. d min 

Remarkably, e\ does not depend on the distance to the 
source if ds 3> dmin- 

We now examine traditional astrophysical scenarios in 
which lensing is or may be observed, and assess whether 
braneworld effects are likely to be detectable. We shall 



see that all of the scenarios have eg <C 1, so the condition 
© is satisfied and our use of the Garriga-Tanaka metric 
is valid. 

Primoridal braneworld black holes may grow to super- 
massive scales by accreting dark energy |28|. Therefore 
let us consider lensing by the supermassive black hole at 
the center of our Galaxy. The black hole has a mass 
of M. = (3.6 ± 0.2) x 10 6 M Q j^, corresponding to 
m. = 5.3 x 10 11 cm = 1.7 x 10~ 7 pc, and a distance 
d.L = 7.9 ± 0.4 kpc [3(J. Assuming that the lens-source 
distance dLS is much smaller than the observer-lens dis- 
tance dL, eq. (|o31) yields 



ej w 6.1 x 10~ 35 



I 



0.2 mm J V pc / 



(66) 



It would be very challenging to measure braneworld ef- 
fects in lensing by the Galactic black hole. 

To increase e\ we need to lower the black hole mass. 
Considering stellar-mass black holes brings us into the 
realm of microlensing (e.g., [ia |3l], |32( ) ■ In Galactic mi- 
crolensing the typical distances between observer, lens, 
and source are kiloparsecs, so eq. (|64|l yields 



2.2 x 10" 31 x 



0.2 mm 



Mm 



dL dLS 



(67) 

where the distances di are in kpc. Compared with the 
example of the Galactic supermassive black hole, we have 
gained six orders of magnitude in mass but lost three in 
distance, so the braneworld effects are still very small. 

An intriguing microlensing system is the binary pul- 
sar J0737— 3039 |3^|. The binary orbit is seen nearly 
edge-on, so that when one neutron star passes behind 
the other in projection there may be significant lensing 
effects |34j. Although the lens in this case is not a black 
hole, it is still a compact object to which our formal- 
ism applies. The Einstein radius in this system is 2550 
km while the gravitational radius is m. = 1.85 km [34j. 
so the weak-deflection regime applies. The lens-source 
distance is given by the semimajor axis of the orbit, 
dLs ~ a = 8.79 x 10 5 km, and is much smaller than 
the distance to the lens, so eq. (|65|l yields 



e} w 6.1 x 10~ 21 x 



0.2 mm 



(68) 



Braneworld effects are still negligible. 

We conclude that in traditional lensing scenarios with 
black holes that are stellar mass or larger, the lensing 
scale is simply much too large to provide an effective 
probe of braneworld effects on the scale I < 0.2 mm. 



VII. ATTOLENSING BY PRIMORDIAL 
BRANEWORLD BLACK HOLES 

To obtain larger braneworld effects, we need to con- 
sider even smaller black hole masses. In the braneworld 
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model, black holes may be created in the early universe 
and survive to the present day with masses as small as 
~ 1kg ~ l(r 30 M Q [H Primordial braneworld 

black holes could be detectable via lensing in the wave 
optics limit, in a phenomenon we call attolensing. 

Attolensing is descended from the concept of fem- 
tolensing in general relativity, introduced by Gould |35j | . 
A black hole of mass M, ~ 10~ 16 M placed at a cosmo- 
logical distance creates images with an angular spacing 
of order a femto-arcsecond. Femtolensing could produce 
observable interference fringes in the energy spectrum of 
a gamma-ray burst at energies in the range of keV to 
MeV. Various aspects of the interference patterns pro- 
duced by f emtolen sing have been studied for general rel- 
ativity [EES Sill 

Likewise, attolensing by primordial braneworld black 
holes would be observable through wave optics effects. 
There are, however, some notable differences between 
femtolensing in general relativity and attolensing in 
braneworld gravity. First, at a given black hole mass, 
braneworld gravity and general relativity predict slightly 
different interference patterns (see Section 0. Second, 
braneworld gravity allows lower mass black holes to sur- 
vive to the present day (compared with GR); as a result, 
there are energy scales at which interference effects would 
be observable for braneworld gravity but not general rel- 
ativity. Third, as a point of terminology, Gould's term 
referred to the scale of the angular image separation. By 
contrast, we choose a term that denotes the mass scale, 
which for wave optics is much more fundamental than 
the angular image separation. As noted above, primor- 
dial braneworld black holes can have an enormous range 
of masses, down to ~ 10 -30 M Q . In quantitative exam- 
ples we take M ~ 10~ 18 A/q to be illustrative (hence 
"attolensing"), but always quote the mass scaling. 

The gravitational radius and time scale for attolensing 
are 



m. = 1.5 x 10 



-13 



M. 



t e = 2.0 x 10~ 23 x (l + z L ) 



M. 



(69) 
s, (70) 



where to be general we consider that the lens may have 
cosmological redshift zl. The bright interference fringes 
appear at energies (see eq. E 



Ef = 210xl±l/l— i_ 
3 l + z L A + Ai 



M. 



1 - 



A 



A + Ai 



MeV 



(71) 



The dark fringes are found by replacing j + 1/4 with 
j + 3/4 (cf. ea. 161 jl. The dimensionless factors involving 
Ao and Ai depend on the source position as shown in 
Figure [5] Braneworld effects introduce a fractional shift 
in the fringe energies (relative to GR) that is given by e\ 
times a factor that is approximately 0.5 for all relevant 
source positions. 




FIG. 2: Dimensionless factors in the energies of attolensing 
interference fringes (see eg. 1711 . which depend on the scaled 
source position ft. 



A. Braneworld Black Holes in the Solar System 

Since the amplitude of braneworld effects decreases 
with the distance to the lens (see eg. 16511 . we can estimate 
the maximum realistic effects by considering primordial 
braneworld black holes that are as close as possible to 
Earth. To estimate how close this might be, let pbh be 
the mean mass density in primordial braneworld black 
holes in the Solar neighborhood. We might take this to 
be some fraction /bh of the density pdm of dark mat- 
ter in our region of the Galaxy. If all braneworld black 
holes have the same mass AT,, then their number density 
is ribh = Pbh/A/,, so the typical distance between black 
holes is 



-1/3 



'bh 



( fbhPdr, 



-1/3 



(72) 



Detailed modeling of our Galaxy indicate that the den- 
sity of dark matter in the Solar neighborhood is pdm = 
(0.011 ± 0.005) M Q pc" 3 HJ. This dark matter density 
yields 



i L ~ 0.93 x / bh 1/3 



M. 



10~ 18 M, 







1/3 



AU. 



(73) 



Note that the natural unit here is Astronomical Units. 
In other words, if M,/M Q < 10~ 13 /bh then the nearest 
primordial black holes reside within our Solar System! 

Put another way, a simple estimate of the total mass 
in primordial black holes within the Solar System is 



bh 



'-nR'p Pbh ~ 3.3 x 10~ 13 / bh M, 



© i 



(74) 
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where Rp w 40 AU is the radius of Pluto's orbit. This 
corresponds to a total number 



bh 



3.3 x 10 5 x /, 



bh 



M 



10- 18 M, 



(75) 



of primordial braneworld black holes in the volume inte- 
rior to Pluto's orbit. 

The angular Einstein radius and braneworld correction 
scale for a primordial braneworld black hole in the Solar 
System are 



§ E ~ 4.1 x 10~ 8 x 



M. 



1/2 



1O- 18 M 
e? - 4.5 x 10~ 5 x 



AU 
M. 



-1/2 



arcsec, 



AU 



(76) 



(77) 



In other words, braneworld corrections for lensing by pri- 
mordial braneworld black holes in the Solar System are 
small but not absurdly so. 

The primordial braneworld black holes we are consid- 
ering are much less massive than asteroids. In fact, the 
total mass in (I74|l is smaller than many large asteroids 
|4l| . and is well within upper limits on dark matter in 
the Solar System derived from Solar System dynamics 
|42L 143] ] . Our estimate of the mass and number of pri- 
mordial black holes in the Solar System is based on the 
assumption that the dark matter is distributed uniformly 
in the Solar neighborhood, which seems plausible because 
stars and the star formation process are inefficient at cap- 
turing dark matter [4^ . 



B. Braneworld Black Holes in the Cosmos 

Rather than considering braneworld black holes in our 
Solar System, we may imagine them spread through- 
out the universe. The angular Einstein radius and 
braneworld correction scale for a primordial braneworld 
black hole at a cosmological distance are 



2.9 x 10~ 15 x 



M. 



2.2 x 10~ 19 x 
2 



1/2 / . \ 1/2 



cLl ds 
M, 



(78) 



(79) 



0.2 mm/ V10- 18 M ( 



'3 



dh d LS 



where the distances are all in Gpc. In other words, both 
the angular image separation and the braneworld cor- 
rections would be difficult to measure for cosmological 
primordial braneworld black holes. However, the inter- 
ference fringes are still very similar to the Solar System 
case; they shift only by the order unity factor 1 + zl (see 
eq.ETJ. 



An obvious question is whether the probability of cos- 
mological attolensing is high enough to be interesting. 
Suppose that primordial braneworld black holes con- 
tribute a fraction fi. of the total energy density of the 
universe. The lensing optical depth — defined to be the 
fraction of the sky covered by Einstein radii, which is 
useful not only for traditional lensing but for attolens- 
ing as well — would then be r such that r/fi» > 0.1, 
with the precise coefficient determined by the distribu- 
tion of source redshifts [3^ Ej] , The optical depth is 
independent of the distribution of primordial black hole 
masses, and is nearly equal to the attolensing probability 
[45| . In other words, if primordial black holes contribute 
a cosmologically significant fraction of matter, then the 
attolensing probability is not very small. 



C. On Observing Attolensing 

For our example black hole of mass M, — 10 _18 Mq, 
eq. I|71|) indicates that semi-classical wave optics ef- 
fects would be seen at energies of tens to hundreds of 
MeV. This energy range will soon be accessible with the 
GLAST satellite, scheduled for launch in 2007 0. (At 
present, energies up to 8 MeV can be observed with the 
INTEGRAL satellite H3- There may be important at- 
tolensing effects that could be seen at these lower ener- 
gies, but they would involve physical optics effects and 
it is not yet possible to examine braneworld black hole 
lensing in the physical optics regime [see Section 0.) 

It is instructive to assess whether GLAST could detect 
the fringe pattern in an attolensed gamma-ray burst. At 
high energies, a typical gamma-ray burst has a power 
law spectrum such that the number of photons per unit 
energy is 



dNo 
dE 



oc E~ 



(80) 



with r\ ss 2.25 (48|. If the gamma-ray burst is attolensed, 
its observed energy spectrum equals the intrinsic spec- 
trum multiplied by the energy-dependent lensing magni- 
fication M. from eq. (|57|l . 



dN, 



obs 



dE 



M(E) E~ 



(81) 



Figure shows a sample energy spectrum for a gamma- 
ray burst attolensed by a braneworld black hole with 
mass M, = 1O~ 18 M0, for a configuration in which the 
scaled source position is j3 = 0.75. The dark fringes 
(valleys) appear to be more prominent than the bright 
fringes (peaks), although this is an artifact of using a 
logarithmic vertical scale. Figure |3Jd then shows an ex- 
ample of a binned spectrum as it might be observed by 
GLAST. In this example, the bin width is about three 
times GLAST's energy resolution, and the errorbars indi- 
cate the statistical uncertainties if there are 100 photons 
detected with energies between 30 MeV and 1 GeV. 
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FIG. 3: Sample energy spectrum of a gamma-ray burst 
attolensed by a braneworld black hole with mass M, = 
10 -18 Mq, for a configuration in which the scaled source po- 
sition is P = 0.75. (a) The dashed lines shows the intrinsic 
energy spectrum, while the dotted line shows the attolensed 
spectrum (both normalized to 100 photons). We follow the 
convention for high-energy spectra and plot the logarithm of 
the number of photons per logarithmic energy interval, (b) 
The dashed and dotted lines reproduce the ideal spectra from 
panel (a). The histogram shows a sample "observation" as 
it might be seen by GLAST. The errorbars indicate the un- 
certainty in the "observed" spectrum if the telescope records 
100 photons with energies between 30 MeV and 1 GeV. (c) 
Cumulative distribution of photon energies. The dashed line 
shows the distribution without lensing; the dotted line shows 
the predicted distribution with lensing; while the solid line 
shows the distribution for a mock observation of 100 photons. 
A Kolmogorov-Smirnov test rejects the hypothesis that this 
gamma-ray burst was not lensed at more than 99% confidence. 



The "observed" spectrum shows an excess near 30-40 
MeV, and a clear dip near 100 MeV, relative to the un- 
lensed spectrum. To determine the significance of such 
features, and more generally to assess the ability to dis- 
tinguish between the unlensed and attolensed cases, we 
use the Kolmogorov-Smirnov test to compare the "ob- 
served" spectrum to the unlensed power law spectrum. 
(The KS test applies to the cumulative distribution of 
photon energies, shown in Figure |3J:.) For this combina- 
tion of mass and source position, 100 counts is typically 
sufficient to reject the hypothesis that the gamma-ray 
burst was not lensed at more than 99% confidence. While 
the number of counts required to reach 99% confidence 
varies with the black hole mass, and to a lesser extent 
the source position, it is important to see that for a typi- 
cal scenario strong results could be obtained with as few 
as 100 counts. GLAST is expected to see of order 20 
gamma-ray bursts per year with more than 100 counts 
above 30 MeV E3. 



D. On Identifying Primordial Braneworld Black 
Holes 

Explicit braneworld effects (characterized by ef) are 
vanishingly small for cosmological scenarios, although 
they might be measurable in Solar System scenarios. 
This does not limit our ability to test braneworld grav- 
ity, though. The reason is that the braneworld model 
makes very different predictions, compared with general 
relativity, about the formation and survival of primor- 
dial black holes |13l Il4 , Il5| . The starkest difference 
is that in GR all primordial black holes smaller than 
~ 10~ 19 M Q would have evaporated by the present day 
(l2^. while in braneworld gravity black holes as small as 
~ 1 kg ~ 10~ 30 Mq may be able to survive to today 
[Tolll6| ]. Therefore, a clear detection of a black hole with 
mass < 10 _19 Mq would violate the GR prediction and 
support the braneworld model. 

Observing attolensing fringes would not only reveal a 
primordial black hole but also place an important up- 
per limit on its mass. From eq. the energies of 
the bright and dark fringes are given by (1 + zl) 
times a dimensionless factor of order unity that depends 
on the source position, (3. The source position could 
be determined from the fringe amplitudes (see eq. I63fl . 
Thus, analyzing the fringes would fix the combination 
M. = (1 + zl)M,. While it may be difficult or impos- 
sible to determine zl, we would know that zl > and 
hence M. < M m . Since we would be seeking evidence 
that there are black holes below some mass threshold, 
upper limits available from attolensing would be useful 
and important. 



VIII. CONCLUSIONS 

We have presented a rigorous and comprehensive an- 
alytical formalism for gravitational lensing due to a 
braneworld black hole described by the Garriga-Tanaka 
metric. Using invariant quantities, we calculated all the 
fundamental geometric optics lensing observables. We 
then used these results to pursue a new direction in 
braneworld black hole lensing: wave optics. We com- 
puted the total magnification in the limit of semi-classical 
wave optics, and gave explicit formulas for the locations 
of the bright and dark interference fringes in the energy 
spectrum of a source lensed by a braneworld black hole. 

Applying our results to realistic examples of lens- 
ing indicates that traditional lensing scenarios involv- 
ing stellar-mass or supermassive black holes will be un- 
able to test braneworld gravity in the foreseeable fu- 
ture. However, attolensing by primordial braneworld 
black holes does provide a powerful opportunity to probe 
braneworld effects via interference fringes in the energy 
spectra of gamma-ray bursts. If primordial braneworld 
black holes contribute a non-negligible fraction of the 
dark matter, there should be many within our Solar Sys- 
tem. These nearby primordial black holes could be used 
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to test braneworld gravity directly by looking for the e\ 
correction terms in the interference pattern, and indi- 
rectly by looking for primordial black holes smaller than 
the evaporation limit predicted for general relativity. Pri- 
mordial black holes would also be spread throughout the 
cosmos. They may produce a measurable probability for 
attolensing of gamma-ray bursts, which would again af- 
ford the possibility of detecting black holes smaller than 
the GR limit. 

It is worth reiterating that a population of primordial 
braneworld black holes in the Solar System is not ruled 
out by current dynamical constraints p2l |43[ . In the 
future, the Laser Interferometry Space Antenna (LISA) 
may be able to detect gravitational impulses from pass- 
ing primordial black holes, but LISA's sensitivity and 
noise from the Moon limit the detectability to M, > 
1O _19 M |5(|. Cosmologically, the lack of obvious fem- 
tolensing in a sample of 118 gamma-ray bursts places a 
weak upper limit on the cosmological density of primor- 
dial black holes in the mass range 10~ 16 -10~ 13 M |5lj : 
but the current constraint neither excludes fractions as 
high as O. ~ 0.1 nor probes masses smaller than about 
10 -16 Mq. In summary, current gravitational constraints 



do not rule out a substantial population of primordial 
black holes; for the future, attolensing may be the only 
way to probe the important mass scale below 10~ 19 Mq. 

While our analysis has been as rigorous and compre- 
hensive as possible, there is one significant limitation. 
Since the full metric describing the spacetime around a 
braneworld black hole is still unknown, we cannot extend 
our analysis of lensing to the physical optics regime at the 
present time. We hope that this new fundamental test 
of braneworld gravity will motivate further attempts to 
determine the full metric for a braneworld black hole. 
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